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We investigate the properties of many-electron systems in two-dimensional polygonal (triangle, 
square, pentagon, hexagon) potential wells by using the density functional theory. The development 
of the ground state electronic structure as a function of the dot size is of particular interest. First 
we show that in the case of two electrons, the Wigner molecule formation agrees with the previous 
exact diagonalization studies. Then we present in detail how the spin symmetry breaks in polygonal 
geometries as the spin density functional theory is applied. In several cases with more than two 
electrons, we find a transition to the crystallized state, yielding coincidence with the number of 
density maxima and the electron number. We show that this transition density, which agrees 
reasonably well with previous estimations, is rather insensitive to both the shape of the dot and the 
electron number. 
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I. INTRODUCTION 

The research of nanoscale electronic structures has 
been expanding continuously. Quantum dotsa represent 
basic electron systems that have been fabricated using 
semiconductor materials for almost fifteen years. Be- 
cause the confinement of electrons in quantum dots or 
'artificial atoms' can be varied at will, they have became 
a playground in which the basic physics of interacting 
electrons can be surveyed and theoretical models can be 
tested. 

In quantum dots the correlation effects between elec- 
trons have to be considered carefully because the exter- 
nal confinement is remarkably weaker than in real atoms, 
where the independent electron model with mean-field 
theories usually gives good results. As the confinement 
strength is lowered, the mutual Coulomb interaction be- 
comes gradually dominant and at a certain point, the 
electron density begins to exhibit localization to classi- 
cal positions in order to minimize the interaction. This 
phenomenon corresponds to the Wigner crystallization 
(WC) in a two-dimensional electron gas (2DEG) (Ref. 
2). According to quantum Monte Carlo simulations, the 
crystallization occurs when the 2D electron density n has 
decreased such that the density parameter r s > 37. Here 
n = l/(%r s ) and r s is given in units of the effective Bohr 
radius a* B — h 2 e/m*e 2 , where e is the dielectric constant 
and m* is the effective electron mass, specific to the semi- 
conductor material in question. In 2DEG with impurities 
the broken translational invariance has been shown to re- 
sult in the WC at a much larger electron density with 
r s ~ 7.5 (Ref. 3). In quantum dots the transition to the 
WC hasJpeen predicted to occur at even higher electron 
densitiesoB. One of the questions in this context is how 
the shape and the electron number of a two-dimensional 
quantum dot affect the crystallization. 

The weak-confinement limit in quantum dots has been 
studied with various. theoretical methods, including-ex- 
act diagonalizationoEl, quantum MonterC-arlo (QMC)QO, 
and unrestricted Hartree-Fock (UHF)ErEI techniques, as 



well as the spin density functional theory (SDFTjffll!. 
In this regime the SDFT allows the formation of spin 
density waves (SDW), i.e. the breaking of the spin 
symmetry^, leading to a lower total energy in the sys- 
tem. In our earlier work for a parabolic six-electron 
quantum dot, we examined the energy difference between 
the polarized and paramagnetic spin states, and showed 
that the SDW solution agrees well with the QMC results, 
in contrast to the symmetry-preserved DFT solutionis. 
However, the problem in the SDFT calculations is the 
fact that only the z-component of the total spin can be 
specified. Therefore, one may ask if a mixed state of 
several eigenstates, corresponding to different S with a 
fixed S z , is physically meaningful, an argument presented 
by Hirose and WingreerO. In fact, a mixed-symmetry 
state is not an eigenstate of the Hamiltonian, but the 
lowest state of a well-defined mixture of symmetries is 
a functionaLof the density at the time of preparation 
of the statecB. The symmetry-broken electron structure 
thus gives more accurate approximations for the energy 
and describes the internal nature of the many-body wave 
function better than the symmetry-preserved solutionEIl 
The criterion for the WC in quantum dots n may be 
determined with several attributes. Egger et alB consid- 
ered three criteria yielding similar results in their QMC 
analysis for parabolic quantum dots. They observed the 
electron density in real space and searched for the con- 
finement at which the shell structure began to form. In 
addition, they monitored a quantity depending on the 
pair-correlation function, and changes in the energy spec- 
tra. Localization may also be observed directly by exam- 
ining|-probability densities of single electronslij. Crefneld 
et aln have studied the systems of two electrons confined 
by polygonal 2D infinite-barrier wells, and their criterion 
for the onset of the WC is the appearence of a local den- 
sity minimum at the center of the dot. In the (S)DFT 
calculations the criterion should be based on the density 
(spin densities). However, it became evident in our study 
that the criterion by Crefneld et al. can not be applied 
for polygonal dots containing more than two electrons, 
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because the electron density has maxima at the corners 
of the dot also at very high electron densities. 

Quantum dots are usually modeled by restricting a cer- 
tain number of electrons to a 2D plane and assuming the 
confining potential to have a parabolic shape. An ex- 
ample of more general modeling is the above-mentioned 
exact diagonalization-atudy by CrefReld et ala. More- 
over, Akbar and Leell3 used the SDFT to study square 
quantum dots which have a small finite extent in the 
third dimension perpendicular to the square. In the case 
of two electrons, they found a good agreement with the 
results by Creffield et al. For two- and four-electron dots, 
Akbar and Lee estimated the onset of the WC to occur 
at r s ~ 6. 

In the present work we employ the SDFT to inves- 
tigate the properties of two-dimensional quantum dots 
with a general polygonal confinement. We concentrate 
on the WC in the weak-confinement limit, which is ob- 
tained simply by increasing the side length of the dot. 
In the numerical calculations we apply a recently devel- 
oped real-space approacbO. As a symmetry-unrestricted 
method it is flexible regarding the applied geometry and 
allows also SDW solutions. In the regime of the spin 
symmetry-broken solutions, we find that further weaken- 
ing of the confinement leads to electron densities with as 
many maxima as there are electrons in the system. We 
show that the appearance of this behavior can be used, 
at least in several cases, consistently as a criterion for the 
onset of the WC for quantum dots of various shapes and 
different electron numbers. 

The outline of the paper is as follows. In Sec. |l| 
we present the theoretical model and the computational 
method of our calculations. From the results we first 
compare the DFT, i.e. spin-compensated calculations for 
a two-electron dot to the exact diagonalization results. 
Then we employ the SDFT and present the symmetry- 
broken solutions first for a two-electron dot, and there- 
after for larger systems. A summary and discussion are 
given m Sec[V[ 



II. METHODS 

The quantum dot material is chosen to be GaAs with 
the effective electron mass m* = 0.067m e and the di- 
electric constant e = 12.4. The effective Bohr radius a* B 
is thus 9.79 nm. The Hamiltonian of a many-electron 
system in a polygonal potential well is written as 



The effective mass approximation (EMA) used with 
the Hamiltonian (Q) represents an alternative to the 
constant-interaction modefia, in which the Coulomb in- 
teraction between the electrons is assumed to be inde- 
pendent of the electron number N. The EMA has been 
shown to be a reliable approximation if the confinement 
is not particularly stronger. 

In the SDFT formalism, the electron density is,*«lved 
self-consistently with the Kohn-Sham equationsE3Eil. To 
approximate the exchange-correlation energy functional, 
we use the local spin density approximation (LSDA) with 
the interpolation form by Tanatar and CeperleyE3 for the 
2DEG. Within the EMA, the single-particle Schrodinger 
equation of the Kohn-Sham scheme reads as 



2m 



:V 2 + KffW 



ipi <a (r) = eiip ia (r), (3) 



where V° s is the effective potential for spin a contain- 
ing the external potential and the Hartree and exchange- 
correlation potentials of the electron-electron interac- 
tions. In the spin-compensated calculations (equal spin 
densities), the scheme reduces to the standard density 
functional theory (DFT) within the local density approx- 
imation (LDA). 

In the self-consistent Kohn-Sham scheme, we perform 
calculations in real space by using two-dimensional point 
grids. The differential equations are discretized with fi- 
nite differences^, and the procedure is efficiently accel- 
erated with multigrid techniquesEa to solve the Poisson 
and the single-particle Schrodinger equations. Applying 
the multigrid scheme in the latter case is a fairly com- 
plicated task because both the eigenfunctions and the 
eigenvalues have to be solved simultaneously. In order 
to avoid nonlinearity problems, the Rayleigh Quotient 
Multigrid (RQMG) methodB is used for the solution of 
the eigenpair corresponding to the lowest eigenvalue. We 
employ this method with a recent generalization to an 
arbitrary number of lowest eigenenergy statesEZI. The dis- 
cretized eigenvalue equation is solved by minimizing the 
Rayleigh quotient (ip\H \tp) / ^1^) on the finest grid, us- 
ing the coarser grids to remove the lower frequency com- 
ponents of the error. The technique reduces remarkably 
the number of relaxation sweeps needed for solving the 
Schrodinger equation. Other advantages of the real-space 
solver are its flexibility with the boundary conditions and 
good suitability for parallel computing. 
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where the external potential has a simple form 

f 0, in the dot 

V ext {x,y) = < 

oo, elsewhere. 



III. WIGNER CRYSTALLIZATION OF TWO 

(1) ELECTRONS 

First we perform DFT calculations on two-electron 
polygonal quantum dots by setting the spin densities 
equal (n-f = n\) to prevent the breaking of the spin 

(2) symmetry. The ensuing total energy with its decom- 
position to Coulomb, kinetic and exchange-correlation 
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energies are given in Fig. [j] for a square dot. As pre- 
dicted, the Coulomb energy becomes increasingly more 
dominant over the kinetic energy as the side length L of 
the dot is enlarged and the WC is expected to occur. 




1 ^0 100 150 200 250 300 350 400 
Side length [nm] 

FIG. 1. Energy composition in a square two-electron quan- 
tum dot as a function of the dot size. 

The ground state electron density distributions for the 
triangular, square, pentagonal, and hexagonal dots at 
three side lengths L = 50, 100, and 400 nm are presented 
in Fig. ||. In the small dots the electron density is lumped 
at the center, whereas the large dots represent Wigner- 
molecule-like behavior, the density being localized near 
the corners in order to minimize the dominating Coulomb 
interaction. The localization is seen to. depend strongly 
on the area of the dot. Creffield et aln defined the sys- 
tem to be a Wigner molecule, when a local minimum first 
appears at the center. According to our calculations, this 
occurs, for instance, in the square as the side length is 



about 100 nm (~ 10a* g ), which agrees with the exact di- 
agonalization results. This qualitative consistency estab- 
lishes the applicability of the density functional approach 
to small systems considered in this study. 

We define the density parameter as r s — yj A/Ntt, 
where A is the area of the polygon. In the case of n cor- 
ners and a side length L we thus get r s — ^-\fj^ cat zT. 
By applying the criterion presented by Creffield et ala for 
the WC transition point, we find the value of r s ~ 3 in all 
four geometries. Akbar and LeeE3 employed the SDFT 
for square 2D quantum dots with an additional harmonic 
confinement along the z-axis. They used a more rigorous 
criterion for the WC, i.e. the breaking of connections 
between the density maxima, and estimated the criti- 
cal value of r s ~ 6 for the transition point. In spite of 
the very different estimation technique, this result is in a 
qualitative agreement with our r s , ~ 3. 

Intuitively, the localization of two electrons into all 
the corners of a polygonal potential well might first ait. 
pear as a slightly odd result. Jefferson and HauslerEEl 
have explained the phenomenon with effective charge- 
spin models. They suggested that the low-energy man- 
ifold of a system of strongly correlated electrons can be 
described properly with an extended single-band Hub- 
bard model. For example, in a square two-electron dot 
the tV-Hamiltonian transforms into the following effec- 
tive Hamiltonian, 

Beff = E + (Ae l2 *i?, /2 + h.c.) , (4) 

where rotates the electrons at opposite corners on a 
diagonal by ir/2. The electron pair may thus tunnel be- 
tween the ground state configurations with an amplitude 
modulated by a factor e l2 *. This explains the four-peak 
structure of the electron density in the Wigner limit, pre- 
dicted already by BryantO. Diagonalization of if e ff gives 
a good approximation for the gpaund state energies ob- 
tained from the tV-HamiltonianEa. 
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FIG. 2. Electron densities in polygonal two-electron quantum dots with different sizes. In the square, pentagon, and hexagon 
the amplitudes have been multiplied by a factor 2. 



IV. SYMMETRY-BROKEN SOLUTIONS 



A. Two-electron dot 



Next we perform the same calculations as above but 
without the restriction n-j- = rij, and consider still the 
ground state solution, for which S z = 0. Comparison of 
the new total energies with the spin-compensated results 
as a function of the dot size reveals an interesting transi- 
tion to a lower energy state. At this point, representing 
already a Wigner-crystallized distribution, the spin sym- 
metry breaks and the result is a SDW-like ground state. 

The relative energy differences between the spin sym- 
metric and SDW-like solutions, corresponding to our 
DFT and SDFT calculations respectively, are shown in 
Fig. H for all the considered geometries. In the triangular 
well the transition to the symmetry-broken ground state 
occurs at a remarkably smaller size than in the other 
three geometries. More precisely, for the triangle we get 
the transition at r s ~ 3.5 and for the square, pentagon, 
and hexagon at r s ~ 4.5. 

In order to explain this behavior, one may first ex- 
amine the lowest one-electron energy states, shown in 
Fig. ^ for the triangular and square quantum dots in the 
symmetry-broken SDFT ground state, as well as in the 



symmetry-preserved DFT solution. In the latter state, 
the three-fold geometry produces more low-lying degen- 
erate levels in the triangle than the four-fold geometry in 
the square. In the SDFT calculations these degeneracies 
are split such that the energy levels become pronouncedly 
spread in the triangular geometry, whereby the lowest 
levels are pushed more efficiently downwards in the tri- 
angle than in the square. There is also a qualitative dif- 
ference between the symmetry-broken electron densities 
in these geometries. As shown in Fig. |^, the spin-up and 
spin-down densities are totally separated in the square, 
whereas in the triangle they share a corner. In the trian- 
gular geometry, the breaking of the spin symmetry can 
thus lower the energy via the exchange-correlation and 
Coulomb contributions relatively more, and with a rel- 
atively smaller cost in the kinetic energy than in the 
square. Nevertheless, in none of these geometries the 
breaking of the spin symmetry enlarges the Fermi gap, 
contrary .to the SDW formation in large, parabolic quan- 
tum dotM 
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FIG. 3. Total energy differences between the DFT (£a) 
and SDFT (Eq) solutions in polygonal two-electron quantum 
dots of four geometries. 
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FIG. 4. Lowest energy levels of triangular and square 
two-electron quantum dots at r s ~ 8. 



The composition of the energy difference between the 
symmetry-preserved and the symmetry-broken states are 
presented in Fig. || for a square dot. Naturally, the 
change in the exchange-correlation energy favors and the 
change in the kinetic energy opposes the transition. The 
behavior of the Coulomb energy is interesting: its strong 
decrease actually initiates the breaking of the spin sym- 
metry. However, as the dot is made larger than L ~ 250 
nm, the Coulomb energy is higher in the SDW-like than 
in the symmetry-preserved state. The phenomenon can 
be understood by having a further look at the electron 
density distributions in the square (Fig. ||). In the SDFT 
solution, the electron density is shifted from the region 
between the opposite spin directions towards the corners. 
At small dot sizes this decreases the Coulomb repulsion 
between the charge peaks in the adjacent corners more 
than the repulsion increases inside the peaks. At large 
distances the opposite is true. 
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FIG. 5. Composition of the energy difference between the 
DFT (Ei) and SDFT (Eo) solutions of a square two-electron 
quantum dot as a function of the dot size. 
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FIG. 6. Difference in the electron densities between the 
DFT (up) and SDFT (down) solutions in a square and tri- 
angular two-electron quantum dot at L — 400 nm. The spin 
alignments are shown in the SDFT case. 
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B. TV > 2 

Then we consider some special cases with more than 
two electrons. The next geometry-independent magic 
configuration after TV = 2 is a six-electron dot. It rep- 
resents an interesting point of comparison with the re- 
sults obtained for a parabolic quantum dot in the weak- 
confinement limit. Wc find that the spin symmetry is 
broken at r s ~ 3.8, 3.1, 4.6, 4.9 in the triangle, square, 
pentagon, and hexagon quantum dots, respectively. 

In the parabolic dot with V ex t(r) = ^uo^r 2 , the r s pa- 
rameter can be estimated from lo\ = e 2 /(eneoem*r^^/N) 
(Ref. 12). Jh pursuance of our earlier work for this 
quantum dotlij, the SDW formation was not found until 
r s ~ 6.6. The sharp corners in the confinement seem 
thereby favor the transition to the symmetry-broken 
state. In the six-electron case, however, the triangular 
geometry is more stable against the transition than the 
square one. A square with TV — 6 represents an incon- 
venient combination, similar to the triangle with TV = 2, 
in which the electrons cannot be evenly divided to the 
corners of the polygon. As the number of the corners 
increases further, the transition shifts to higher r s val- 
ues, approaching the point of the SDW formation in the 
parabolic quantum dot with a circular symmetry. 

For TV = 6, we consider also the possibility of spin- 
polarization, i.e. the S z — 3 state becoming the ground 
state in the low density limit. The energy differences be- 
tween the polarized (S z = 3) and paramagnetic (S z = 0) 
states for triangular and square geometries as a function 
of r s are shown in Fig. ^. For comparison, the SDFT re- 
sults for the parabolic quantum dottil are also presented, 
the latter showing spin-polarization at r s > 12. We were 
not able to obtain well-converged results for the triangle 
and square quantum dots at large r s values. Therefore 
we can only speculate by extrapolation that the polariza- 
tion could occur in the triangle and square slightly earlier 
than in the parabolic quantum dot. 

Besides the geometry, we can study how the number of 
electrons affects the breaking of spin symmetry. First we 
consider a square dot with TV = 6, 8, and 12, which all 
correspond to completely filled shells. Fig. || shows the 
energy difference between the spin symmetry-preserved 
and -broken solutions as a function of the r s . For TV = 6 
and TV = 12, the spin symmetry breaks at r s ~ 1.7 and 
r s ~ 1.1, respectively, whereas the ground state of the 
TV = 8 dot remains spin symmetric until r s ~ 2.8. How- 
ever, the energy difference grows rapidly in this dot, being 
considerably larger than in the TV = 6 dot at r s ~ 10. 
In the large dots the SDFT solutions show pronounced 
localization of the spin densities as can be seen in Fig. 
|9|. For TV = 6 and 8, the number of maxima in the to- 
tal electron densities equals to the number of electrons, 
leading to tt and 7r/2 rotational symmetries in these sys- 
tems, respectively. The spin symmetry can be considered 
to be broken more completely in the TV = 8 dot, where 
the density peaks with the same spin are located on di- 



agonally opposite vertices, in contrast to the TV = 6 dot 
where they lie on adjacent corners. The interaction is 
thus minimized more efficiently in the TV = 8 dot, corre- 
sponding to a relatively rapid decrease of the total energy 
shown in Fig. ^|. For comparison, the results for a tri- 
angular quantum dot with TV = 6 are also presented. In 
that system, the increase in the energy difference resem- 
bles the behavior of the TV = 8 square dot, reflecting a 
similar symmetry-broken geometry (see Fig. O below). 
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FIG. 7. Total energy differences between the S z = 3 and 
S z = states in the triangular (triangle markers), square 
(square markers), and parabolic (dashed line) six-electron 
quantum dots. 
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FIG. 8. Total energy differences between the DFT {Ex) 
and SDFT (Eo) solutions in square quantum dots with TV = 
6, 8, and 12 (square markers). Results for the TV = 6 triangle 
quantum dot are also given (triangle markers). 
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FIG. 9. Electron densities of the DFT (up) and SDFT 
(down) solutions in the N = 6 and N = 8 square quantum 
dots with side lengths L = 300 nm. The spin alignments are 
shown in the SDFT case. 



After the breaking of the spin symmetry, there can 
be seen only four density maxima in the corners of the 
N = 6 and N = 8 square quantum dots, resembling 
the DFT solution (the upper row of Fig. ||). The dot 
size has to be increased substantially before the maxima 
in the middle of the edges appear (the lower row). We 
can observe the same behavior in the N — 6 triangle 
and N = 10 pentagon, in both which the spin symmetry 
breaks at r s ~ 2.3. Their density distributions at large 
r s values are given in Fig. [l0| In all these four cases, 
the number of density maxima equals to the number of 
electrons in the system. Therefore the appearance of the 
last density peaks can be considered as the final stage in 
the onset of the WC in the SDFT formalism. 
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FIG. 10. Electron densities at r a ~ 8 in triangular and pen- 
tagonal quantum dots with N = 6 and N = 10, respectively. 



become closer to each other and are remarkably lowered 
in comparison with the symmetry-preserved DFT solu- 
tion. This condensation occurs similarly in all the dots, 
in which the electron density localizes to a number of 
maxima coinciding with the number of electrons. The 
appearance of the last density maxima thus drives the 
lowest energy levels towards degeneracy. The complete 
degeneracy would be the ultimate state for the Wigner 
crystal. In Fig. [ll] one can also notice that the Fermi 
gap is considerably larger in the symmetry-broken solu- 
tion than in the symmetry-preserved case, resembling the 
situation in large, parabolic quantum dota!3. 
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FIG. 11. Lowest single-particle energy levels in a N = 6 
square quantum dot with side lengths L = 100 and 400 nm. 
Solid and dashed lines correspond to the occupied and un- 
occupied staes, respectively. The levels are nondegenerate, 
except the doubly degenerate levels denoted by the numbers 
(2). 

We have carefully determined the r s values at which 
the last maxima appear and found astonishingly simi- 
lar values for the different systems studied, although the 
breaking of the spin symmetry occurs on a broad r s scale. 
The critical values of r s ~ are 3.8 and 4.0 for N = 6 and 
N = 8 square dots, respectively, and r s ~ 3.9 for both 
the N = 6 triangle and the N = 10 pentagon. In the 
case of two-electron dots, the above criterion for the WC 
cannot be applied, but the onset of the spin symmetry- 
broken state gives an reasonable estimate of r s ~ 3.5 for 
the triangular N = 2 dot and r s ~ 4.5 for the other 
N = 2 polygonal quantum dots. Our estimate of r s ~ 
4.0 for the WC transition point is consistent with the 
results for small, parabolic quantum dota3~113. It is also 
clearly smaller than r s ~ 7.5 obtained fot-thc fluid-solid 
transition in 2DEG containing impuritiesu. 



In order to analyze the appearance of the last density 
peaks, we show in Fig. |ll| the lowest one-electron en- 
ergy levels for the N — 6 square quantum dot with side 
lengths L = 100 and 400 nm, corresponding to r s ~ 2 
and 9, respectively. At the smaller size, the spin sym- 
metry has already broken and split the DFT degenera- 
cies. As the dot is made larger, the two lowest states 



V. SUMMARY 

We have studied the electronic properties of polygo- 
nal two-dimensional quantum dots by employing the spin 
density functional theory. The numerical calculations 
are performed with a symmetry-unrestricted real-space 
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scheme. Especially, we have been focused on the be- 
havior of these systems at the weak-confinement limit, 
where the role of the electron-electron interactions be- 
comes dominating and eventually leads to the formation 
of the so-called Wigner molecules. 

First we have shown that the density functional the- 
ory is capable to reproduce, in agreement with the ex- 
act diagonalization studies, the behavior of the electron 
density in polygonal two-electron quantum dots as the 
spatial size of the potential well increases. 

The spin density functional theory leads inevitably to 
the breaking of the spin symmetry. For different geome- 
tries and different electron numbers, this occurs in a wide 
range of average electron densities or r s parameters. The 
spin symmetry-broken density shows for certain geome- 
tries and electron numbers a gradual transition, such that 
the number of density maxima coincides with the num- 
ber of electrons. We use the the appearance of the last 
density maxima as the criterion for the Wigner crystal- 
lization and obtain r s ~ 4.0 for the critical density. This 
value does not depend strongly on the geometry, nor the 
electron number of the quantum dot, and is in agreement 
with Quantum Monte Carlo results. 
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